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Abstract 

The higher twist corrections h N (x)/Q 2 to the spin dependent proton and 
neutron structure functions (x, Q 2 ) are extracted in a model independent way 
from experimental data on and found to be non-negligible. It is shown that 
the NLO QCD polarized parton densities determined from the data on g\, in- 
cluding higher twist effects, are in good agreement with those found earlier from 
our analysis of the data on g\jF\ and A\ where higher twist effects are negligible. 
On the contrary, the LO QCD polarized parton densities obtained from the data 
on g\, including higher twist, differ significantly from our previous results. 

PACS:13m.Bb; 13.88+e; 12.38.-t;13.30.-a 



1 Introduction 



Spurred on by the famous EMC experiment [1] at CERN in 1987, there has been 
a huge growth of interest in polarized DIS experiments which yield more refined in- 
formation about the partonic structure of the nucleon, i.e., how the nucleon spin is 
divided up among its constituents, quarks and gluons. Many experiments have been 
carried out at SLAC, CERN and DESY to measure the longitudinal (A\\) and trans- 
verse (A±) asymmetries and to extract from them the photon- nucleon asymmetries 
Ai(x,Q 2 ) and A 2 (x, Q 2 ) as well as the nucleon spin-dependent structure functions 
gi(x,Q 2 ) and g 2 (x,Q 2 ). Many theoretical analyses of the world data on A 1 and g x 
based on leading (LO) and next-to-leading order (NLO) calculations in perturbative 
QCD have been performed in order to test the spin properties of QCD and extract 
from the data the polarized parton densities. t It was demonstrated that the polarized 
DIS data are in an excellent agreement with the pQCD predictions for A^(x, Q 2 ) and 
g^(x,Q 2 ). What also follows from these analyses is that the limited kinematic range 
and the precision of the present generation of inclusive DIS experiments are enough to 
determine with a good accuracy only the polarized parton densities (Aw + Au)(x, Q 2 ) 
and (Ad + Ad)(x,Q 2 ). The polarized strange sea density (As + As)(x, Q 2 ) as well 
as the polarized gluon density AG(x, Q 2 ) are still weakly constrained, especially AG. 
The non-strange polarized sea-quark densities Au and Ad cannot be determined, even 
in principle, from the inclusive DIS experiments alone without additional assumptions. 

There is, however, an important difference between the kinematic regions of the 
unpolarized and polarized data sets. While in the unpolarized case we can cut the low 
Q 2 and W 2 data in order to eliminate the less known non-perturbative higher twist 
effects, it is impossible to perform such a procedure for the present data on the spin- 
dependent structure functions without loosing too much information. This is especially 
the case for the HERMES, SLAC and Jefferson Lab experiments. So, to extract the 
polarized parton densities from the experimental data on gi(x,Q 2 ) the higher twist 
corrections have to be included in the data fits. Note that the polarized parton densities 
in QCD are related to the leading- twist expression of g±. 

It was shown [3, 4, 5] that to avoid this problem and to determine polarized parton 
densities less sensitive to higher twist effects it is better to analyze data on 
gi/Fx) using for the g 1 and F 1 structure functions their leading twist (LT) expressions. 

'''Note that the theoretical analyses have been mainly concentrated on the Ai(gi) data because the 
measurements of the quantities ^2(32) are much less accurate with the exception of the very recent 
data of El 55 Collaboration at SLAC [2]. Another reason is that the theoretical treatment of gi is 
much more complicated. 
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It is found that if for (<7i)lt an NLO approximation is used, the "effective higher 
twist" corrections to Ai, extracted from the data, are negligible and consistent with 
zero within the errors, which means that the higher twists corrections (HT) to g± 
and Fi approximately cancel in the ratio gi/Fi, or more precisely, ((?i)ht / '(#i)lt ~ 
(-Pi)ht/ (-Pi)lt- 

In this paper we present a detailed study of the higher twist contributions h N (x) / Q 2 
to the nucleon structure function gi(x, Q 2 ). The quantities h N (x) have been extracted 
from the data in a model independent way. The role of higher twists in the determina- 
tion of the polarized parton densities is discussed. 



2 Connection between Theory and Experiment 

The nucleon spin- dependent structure function g^(x,Q 2 ) is a linear combination of 
the asymmetries AW and A 1 ^ (or the related virtual photon-nucleon asymmetries A± 
and A2) measured with the target polarized longitudinally or perpendicular to the 
lepton beam, respectively. The most direct way to confront the QCD predictions to 
the data is a fit to data on the ratio of the structure functions, g^ / F^ . Such data 
have been directly presented by SLAC/E143 and SLAC/E155 experiments [6]. Most 
of the Collaborations, however, have presented data on the asymmetry A± which, in 
practice, are data on AW /D. The photon-nucleon asymmetry A± and the ratio g± / F^ 
are related to the measured longitudinal asymmetry A^ by 

A N 

A N , = (1) 

n N AN 

(l+7 2 )|v = -^ + (7-^, (2) 

where D denotes the photon depolarization factor, r] and 7 are kinematic factors. 77 is 
proportional to 7 and 7 is given by 

7" = ^- (3) 

In (3) Mn is the nucleon mass. It should be noted that in the SLAC and HERMES 
kinematic regions 7 2 cannot be neglected on LHS of (2). 

The magnitude of A% has been measured by SMC, SLAC/E143 and SLAC/E155 
and found to be small [2, 7]. Then to a good approximation its contribution to the 
RHS of Eqs. (1) and (2) can be neglected and A± and gi / Ff can be expressed as 

A N 

A? = -i, (4) 
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(l + 7 2 )|v = (5) 

It is important to note that due to the additional small factor (7 — 77) in (2) the ratio 
g^/F^ is better approximated by the measured asymmetry (Eq. (5)) than the 
virtual photon- nucleon asymmetry (Eq. (4)). 

Using (4) and (5) we reach the well known relation 



\N(„ r\2\ ~ it 1 , 2\ 9ii x -> Q 2 ) 

F ± N (x,Q 2 ) 



Af(x,Q 2 ) - (l + 7 2 )|ye^r (6) 



usually used in the literature. However, as was already mentioned, we have to keep in 



exp 5 



mind that the presented experimental values on A± [1, 8] neglecting A 2 , "A±(x, Q 2 )" 
are really the experimental values of A^(x, Q 2 )/D and that the latter quantity is very 
well approximated by (1 + 7 2 )(?f r / F^ . 

Using the relation between the unpolarized structure function Fi(x, Q 2 ) and the 
usually extracted from unpolarized DIS experiments F 2 (x, Q 2 ) and R(x,Q 2 ) 

2xF? = F 2 N (1 + 1 2 )/(1 + R N ) (N = p,n,d) (7) 

Eq. (6) can be rewritten as 

Q 2 ) = 9 L M 2m i + Q 2 )] ■ (8) 



Up to now, two approaches have been mainly used to extract the polarized parton 
densities (PPD) from the world polarized DIS data. According to the first [4, 5] the 
leading twist LO/NLO QCD expressions for the structure functions g± and Ff have 
been used in (6) in order to confront the data 



9i{x,Q 2 ) 
Fi(x,Q 2 ) 



exp 



gi(x,Q 2 ) LT 
F 1 (x,Q 2 ) LT ' 



1 ,, ,-t i-,9i( x :Q 2 )lt 

A^x.Q)^ & (1 + 7) ^7—7^ • (9) 

In (9) we have dropped the nucleon target label N. It was shown [3, 4, 5] that in this 
case the extracted from the data "effective" HT corrections h Al (x) to A 1 

are negligible and consistent with zero within the errors, h Al (x) ~ (see Fig. 1). What 
follows from this result is that the higher twist corrections to g\ and F\ compensate 
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each other in the ratio gi/Fi and the PPD extracted this way are less sensitive to 
higher twist effects. We stress again that the polarized parton densities in QCD are 
related only to the leading- twist part of g\. 

According to the second approach [9, 10], g\/F\ and A 1 data have been fitted using 
phenomenological parametrizations of the experimental data for F 2 (x, Q 2 ) and R(x, Q 2 ) 

~ gi{x,Q 2 ) 
F 1 (x,Q 2 ) 

A 1 (x,Q 2 ) exp & yi )^^ 2 )LT 2x[l + R(x,Q 2 ) exp ] . (11) 

"2 \%i jexp 

Note that such a procedure is equivalent to a fit to (<7i) exp , but it is in principle better 
than the fit to the g\ data themselves actually presented by the experimental groups. 
The point is that most of the experimental data on g\ have been extracted from the A\ 
and gi/Fi data using the additional assumption that the ratio g\fF\ does not depend 
on Q 2 . Also, different experimental groups have used different parametrizations for F 2 
and R. 

In the analyses [11, 12] a procedure has been used which is somehow a mixture 
between the two methods above, but bearing in mind the sensitivity of the results to 
higher twist effects it is analogous to the second one. In these fits the leading twist 
expression for F 2 instead of its experimental values has been used 

' ' 9\{x,Q 2 ) 
F 1 (x,Qi) 

A 1 (x,Q 2 ) exp & ^l^L 2 x[l + R(x,Q 2 ) exp ] . (12) 

It was shown by GRSV [4] that if the second approach (11) is applied to the data 
(i*2 and R are taken from experiment) allowing at the same time "effective higher 
twist" contribution h Al (x)/Q 2 to the RHS of (11), h Al (x) is found to be sizeable and 
important in the fit. In other words, bearing in mind that a lot of data on A-y and 
<7i/i*i are at small Q 2 special attention must be paid to higher twist corrections to 
the structure function g\. To extract correctly the polarized parton densities from the 
gi data these corrections have to be included into data fits. Note that a QCD fit to 
the data in this case, keeping in gi{x,Q 2 )QCD only the leading-twist expression, leads 
to some "effective" parton densities which involve in themselves the HT effects and 
therefore, are not quite correct. 



9i(x,Q )lt Jl + R(x,Q ) exp] 



( rn 



F 2 (x,Q 2 



-2x- 



I exp 



:i+7 2 ) 



f-'.:r.7) 



gi(x,Q 2 ) hT [1 + R{x,Q 2 ) 

expi 

x + 7^ 
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3 Higher twist effects in gi(x y Q 2 ) 



The usual pQCD expression for the nucleon structure function g{(x,Q 2 ), in terms of 
polarized quark and gluon densities, arises from the contribution of the leading twist 
(r = 2) QCD operators and in NLO has the form (a similar formula holds for g'[ l ): 

£(x, Q 2 ) pQ cd = I £ #(A? + Aq) ® (1 + ^SC q ) + ^AG ® (13) 

where Aq(x, Q 2 ), Aq(x, Q 2 ) and AG(x, Q 2 ) are quark, anti-quark and gluon polarized 
densities in the proton, which evolve in Q 2 according to the spin-dependent NLO 
DGLAP equations. 5C(x) q> G are the NLO spin-dependent Wilson coefficient functions 
and the symbol ® denotes the usual convolution in Bjorken x space. Nf is the number 
of active flavors. In LO QCD the coefficients SC(x) q and 8C(x)g vanish and the 
polarized parton densities in (13) evolve in Q 2 according to the spin-dependent LO 
DGLAP equations. 

It is well known that at NLO and beyond, the parton densities as well as the 
Wilson coefficient functions become dependent on the renormalization (or factorization) 
scheme employed. ^ Two often used schemes are the MS and the JET schemes. Both 
the NLO polarized coefficient functions [13] and the NLO polarized splitting functions 
(anomalous dimensions) [14] needed for the calculation of gi{x,Q 2 ) in the MS scheme 
are well known at present. The corresponding expressions for these quantities in the 
JET scheme can be found in [15]. 

However, there are other contributions to gi, arising from QCD operators of higher 
twist (HT), namely r > 3, which are related to multi-parton correlations in the nucleon. 
It can be shown that these give rise to contributions to gi(x,Q 2 ) that decrease like 
inverse powers of Q 2 . The leading term has the form h N (x,Q 2 )/Q 2 , where h N (x,Q 2 ) 
could have a slow, logarithmic dependence on Q 2 . 

In the kinematic regime where such terms might be relevant it is important for 
consistency to realise that the QCD expression (13) is derived under the assumption 
that Q 2 » Mjq. There will thus be purely kinematic corrections to (13), which involve 
a power series in M^/Q 2 with small coefficients. The leading term of these so-called 
target mass corrections (TMC) therefore has a Q 2 behavior similar to the genuine HT 
terms, but it is not a dynamical HT effect. 

In view of this we shall write 

9i(x,Q 2 ) = gi (x, Q 2 ) LT + gi (x, Q 2 ) HT , (14) 

tOf course, physical quantities such as the virtual photon-nucleon asymmetry A\(x, Q 2 ) and the 
polarized structure function gi(x,Q 2 ) are independent of choice of the factorization convention. 
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where we have dropped the nucleoli target label N. In (14) 

9i(x, Q 2 )lt = m(x, QXqcd + h TMC (x, Q 2 )/Q 2 , (15) 

where h™ c (x,Q 2 ) is exactly calculable [16, 17] and 

gi(x,Q 2 ) UT = h(x,Q 2 )/Q 2 . (16) 

As mentioned above, h(x, Q 2 ) denotes the dynamical higher twist power corrections 
to gi which represent the multi-parton correlations in the target. The latter are non- 
perturbative effects and their calculation is model dependent (see, e.g., [18] and refer- 
ences therein). That is why a model independent extraction of the dynamical higher 
twists h(x) from the experimental data is important not only for a better determination 
of the polarized parton densities but also because it would lead to interesting tests of 
the non-perturbative QCD regime. 



4 Method of Analysis 

In this Section we will briefly describe the method of our analysis of the data on inclusive 
polarized DIS taking into account the higher twist corrections to the nucleon structure 
function gf(x,Q 2 ). In our fit to the data we have used the following expressions for 
g 1 /F 1 and A- L : 

^f(x,g 2 )] 9?(x,Q 2 )ur + h N (x)/Q\ [1 + R(x, Q 2 ) exp ] 



exp - / ex P 



F 2 N (x,Q*) exp ~ (l + 7 2 ) 



Fi(x,Q 2 ) 

A N Ax,Q 2 U P <* ^ fj^ + n ^ (x)/Q2 2x[l + R(x, Q 2 ) exp ] . (17) 

r 2 \X, \q£ )exp 

where gi(x, Q 2 )lt is given by the leading twist expression (15). In (17) h N (x) are 
a measure of the dynamical higher twists. In our analysis their Q 2 dependence is 
neglected. It is small and the accuracy of the present data does not allow to determine 
it. For the unpolarized structure functions F^i^x, Q 2 ) eX p and R(x, Q 2 ) eX p we have used 
the NMC parametrization [19] and the SLAC parametrization R1998 [20], respectively. 

As in our previous analysis [5], for the input LO and NLO polarized parton densities 
at Ql = 1 GeV 2 we have adopted a simple parametrization 

xAu v (x, Q Q ) 

Qo)> 

xAd v (x,Ql) = r] d A d x ad xd v (x,Ql), 
xAs(x,Q 2 ) = r] s A s x as xs(x,Q 2 ), 

xAG(x,Q 2 ) = r) g A g x a °xG(x,Ql), (18) 
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where on RHS of (18) we have used the MRST98 (central gluon) [21] and MRST99 
(central gluon) [22] parametrizations for the LO and NLO(MS) unpolarized densities, 
respectively. The number of active flavors is N f = 3. The normalization factors Aj, 
in (18) are fixed such that rji are the first moments of the polarized densities. To 
fit better the data in LO QCD, an additional factor (1 + ^ v x) on RHS is used for 
the valence quarks. Bearing in mind that the light quark sea densities Am and Ad 
cannot be, in principle, determined from the present inclusive data (in the absence of 
polararized charge current neutrino experiments) we have adopted the convention of a 
flavor symmetric sea 

Au sea = Au = Ad sea = Ad = As = As. (19) 

The first moments of the valence quark densities r\ u and i] d are constrained by the 
sum rules 

a 3 = g A = F + D = 1.2670 ± 0.0035 [23], (20) 

a 8 = 3F - D = 0.585 ± 0.025 [11], (21) 

where 03 and a 8 are non-singlet combinations of the first moments of the polarized 
parton densities corresponding to 3 rd and 8 th components of the axial vector Cabibbo 
current 

a 3 = (Au + Au) (Q 2 ) - (Ad + Ad) (Q 2 ) , (22) 
a s = (Au + Au) (Q 2 ) + (Ad + Ad) (Q 2 ) - 2(As + As) (Q 2 ) . (23) 

The sum rule (20) reflects isospin SU(2) symmetry, whereas (21) is a consequence 
of the SU (3) / flavor symmetry treatment of the hyperon /3-decays. While the isospin 
symmetry is not in doubt, there is some question about the accuracy of assuming 
SU(3)f in analyzing hyperon /3-decays. We have previously studied the sensitivity of 
the polarized parton densities to the deviation of a§ from its SU(3) flavor symmetric 
value (0.58). The results are given in [24]. In this analysis we will use for a 8 its SU(3) 
symmetric value (21). 

In our past papers we have used the Jacobi polynomial method to yield the struc- 
ture functions g^(x, Q 2 )lt from their Mellin moments in n space. The details of this 
procedure are given in [25]. But the accuracy of this method is limited in the low x 
region, x < 0.01, so we have now used the inverse Mellin-transformation method (see, 
e.g., [26]) which reconstructs very precisely g^(x, Q 2 )lt from its moments in the whole 
x region. We have repeated our fits without including HT corrections in order to com- 
pare to our previous results [5] obtained by the Jacobi polynomial method. We have 
found very good agreement between the results obtained by both methods. The reason 
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is that the present kinematic x region of the polarized DIS data coincides with the do- 
main where the Jocobi polynomial method works well. Also, the difference between the 
structure functions calculated by Jacobi and inverse Mellin-transformation methods is 
much smaller than the accuracy achieved in the present polarized DIS experiments. 
Nevertheless, bearing in mind its universality and, in particular, its applicability to the 
semi-inclusive DIS processes, we have decided to use the inverse Mellin-transformation 
approach in this analysis. 

The unknown higher twists h N (x) in (17) have been extracted from the data follow- 
ing the method used in [27] and [28] for the higher twist corrections to the unpolarized 
structure functions F 2 and xF 3 , respectively ^ The measured x region has been split 
into 5 bins and to any rr-bin two parameters and hf 1 ^ have been attached. We have 
found that for a deutron the relation h^f 1 = 0.925(/ij p ' ) + h\ n ^)/2 is a good approxima- 
tion. So, to the parameters connected with the input PPD (18) we add the parameters 
hP and/if , (t = l,2,.,5). 

All free parameters 

{ a u , a d , a„ a g , rj s , 7] g (7„, 7 d ); hf\ } (24) 

have been determined from the best fit to g\jF\ and A\ data using (17), i.e., effectively 
by fitting (gi) e xp- Note that in the calculations of gi(x, Q 2 )lt we have used for the 
strong coupling constant a s (Q 2 ) the same procedure as in our previous paper [5] (see 
the details in LSS2001 FORTRAN code at http://durpdg.dur.ac.uk/HEPDATA/PDF). 

5 Results 

In this section we present the numerical results of our fits to the world data on g\jF\ [6] 
and A\ [1, 8]. The data used (185 experimental points) cover the following kinematic 
region: 

0.005 < x < 0.75, 1 < Q 2 < 58 GeV 2 . (25) 

The total (statistical and systematic) errors are taken into account. The systematic 
errors are added quadratically. 

We prefer to discuss the results of the NLO analysis in the JET (or so-called chirally 

invariant) factorization scheme [30]. In this scheme the first moment of singlet AS(Q 2 ), 

as well as the strange sea polarization (As + As)(Q 2 ), are Q 2 independent quantities. 

'''Note that the moments of the g\ higher twists have been studied in the SLAC/E143 paper [6] as 
well as in [17, 29]. 
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Then, it is meaningful to directly interpret AE as the contribution of the quark spins 
to the nucleon spin and to compare its value obtained from the DIS region with the 
predictions of the different (constituent, chiral, etc.) quark models at low Q 2 . Later we 
will briefly comment on the scheme dependence effects on the results of the analysis. 

5.1 Higher twist effects 

The numerical results of our fits to the data are summarized in Tables 1 and 2. As seen 
from the values of x 2 per degree of freedom Xdf i n Table 1, a very good description of 
the data is achieved. The best LO and NLO(JET) fits correspond to Xdf.lo = 0.892 
and to Xdf,nlo = 0.858. We have found that the fit to the data is significantly im- 
proved, especially in the LO case, when the higher twist corrections to gf are included 
in the analysis (see Table 2). In contrast to the case when the HT corrections to g x 
are not taken into account in the fits, the value of Xlo(HT) is very close to that of 
Xnlo(HT), which is an indication that the tail of the neglected higher order logarithmic 
corrections to g 1 resemble a power behavior of order 0(1/Q 2 ) [31]. A similar behavior 
of x 2 has been observed in the QCD analysis of the unpolarized structure function 
xF 3 {x,Q 2 ) in [28]. 

The extracted higher twist corrections to the proton and neutron spin structure 

functions, h p (x) and h n (x), are shown in Fig. 2. As seen from Fig. 2 the corrections 

for the proton and neutron have a different shape. While h p (x) changes sign in the 

LO as well in the NLO case, h n (x) is non-negative in the measured x region in both 

cases. One can see also that the HT corrections to the proton structure function g\ 

appear to be smaller when for (gi)ur the NLO approximation is used. In Fig. 3 we 

demonstrate how the choice of the factorization scheme for the perturbative calculation 

°f (fl'i)LT influences the higher twists results. The results are presented for the JET 

and MS schemes. It is seen that the HT corrections to g± in both cases coincide within 

the errors. The small deference between the central values could be considered as an 

estimation of the NNLO effects in (gi)ur^- 

Hn [10] the HT terms have been discussed using for them two phenomenological parametrizations. 
The authors conclude that they do not find a significant higher twist contribution to gi in an NLO 
treatment of (<?i)lt- On other hand, studying the higher twist effects in the moments of gi, it was 
shown in [17] that while the first moment of higher twist is quite small, the higher order moments are 
relevant at Q 2 ~ few GeV 2 . 
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5.2 NLO polarized parton densities 

Let us discuss now the polarized parton densities extracted from the data in the pres- 
ence of the HT corrections to g\. We will call this set of parton densities PD(g]" T + HT). 
In Fig. 4 we compare the NLO(JET) polarized PD(g^ LO + HT) with those obtained in 
our analysis [5] where we performed fits to the data according to (9). We will call the 
latter PD(gf LO /F^ LO ). As seen from Fig. 4 the two sets of polarized parton densities 
are very close to each other. This is a good illustration of the fact that a fit to the data 
on gi/F\ ) using for the g\ and F\ structure functions their NLO leading twist 

expressions (xdfnlo = 0.859) is equivalent to a fit to the g 1 data taking into account 
the higher twist corrections to g 1 (Xdf nlo = 0.858). In other words, this analysis con- 
firms once more that the higher twist corrections to g\ and F\ approximately cancel in 
the ratio gi/ F x . 

5.3 LO polarized parton densities 

Let us turn now to the LO polarized parton densities. In LO QCD AG(x, Q 2 ) does not 
contribute directly to g\ and the gluons cannot be determined from DIS data alone. 
For this reason the LO fit to the data was performed using for the input polarized 
gluon density AG(x, Ql) the one extracted in the NLO fit to the data: 

AG(x, Q 2 ) LO = AG(x, Q^)nlo(jet) • (26) 

It is important to note that in the polarized case the LO approximation has some 
peculiarities compared to the unpolarized one. As a consequence of the gluon axial 
anomaly, the difference between NLO anti-quark polarizations Aqi in different factor- 
ization schemes can be quite large, comparable in magnitude to the Aqi themselves 
(see, e.g., [5]). In this case the leading order will be a bad approximation, at least for 
the polarized sea-quark densities. Also, bearing in mind that in polarized DIS most of 
the data points are at low Q 2 , lower than the usual cuts in the analyses of unpolarized 
data (Q 2 > 4 — 5 GeV 2 ), the NLO corrections to all polarized parton densities are 
large in this region and it is better to take them into account. Nevertheless, the LO 
polarized parton densities may be useful for some practical purposes; e.g., for prelim- 
inary estimations of the cross sections in future polarized experiments, etc. They are 
also needed for comparison with those extracted from semi-inclusive DIS data [32], 
where the NLO QCD analysis is still very complicated. The extracted LO polarized 
parton densities PD(g^° + HT) are shown in Fig. 5. Also shown in Fig. 5 are the 
LO polarized PD((yf^°/F 1 LO ) obtained in our analysis [5]. In contrast to the NLO case, 
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the two sets of LO polarized densities are significantly different. As a result we obtain 
different theoretical curves for gi (see Fig. 6). To illustrate how these curves fit the 
data, the SLAC/E143 experimental proton data at Q 2 = 5 GeV 2 are also shown. As 
seen from Fig. 6, the proton structure function g{(x, Q 2 ) LO calculated using the LO 
polarized ~P~D{g\° / F^°) does not agree with the data for x < 0.25. Note that at the 
same time the ratio g\° / F^° fits the world data on g±/Fi and A 1 quite well in the 
measured (x, Q 2 ) region (xdf = 0.921) [5]. The main reason for this peculiarity is that 
the LO approximations for the unpolarized structure functions F 2 (x, Q 2 ) and Fi(x, Q 2 ) 
(or R(x, Q 2 )) presented in the literature are not self consistent. The unpolarized parton 
densities in leading order QCD (more correctly in leading logarithmic approximation 
LLA) are usually extracted including in the data set of the analysis only the experi- 
mental data on F 2 and ignoring the data on R (or Fl). Remember that in LLA of QCD 
the structure functions satisfy the Callan-Gross relation 2xF 1 (x : Q 2 ) = F 2 (x, Q 2 ) which 
leads to R — AM 2 x/Q 2 . However, the experimental data on F 2 and R do not satisfy 
these relations in a large kinematic region. They are approximately satisfied only at 
large x and/or large Q 2 . At small x and Q 2 the experimental values of F 2 are larger 
than those of 2xFx by up to 30%. That is why the extracted sets of LO unpolarized 
parton densities (MRST, CTEQ, etc.) are not quite consistent. While they fit well the 
data on F 2) they badly fail to describe the R{Fi) data in the region of small x and Q 2 . 
One way to improve the situation would be to perform a LO QCD fit including in the 
data set the R{F{) data too. Also, if the data at low Q 2 , lower than the usual cuts 
(Q 2 > 4 — 5 GeV 2 ), are included in the analysis, the higher twist corrections to F 2 and 
Fi should be taken into account. 

6 Conclusion 

We have analyzed the world data on inclusive polarized deep inelastic lepton-nucleon 
scattering in leading and next-to-leading order of QCD including in the analysis the 
higher twist h N (x)/Q 2 and the target mass corrections to the nucleon spin structure 
function g^(x, Q 2 )- We have found that the fit to the data on gi is essentially improved, 
especially in the LO case, when the higher twist terms are included in the analysis. 
The x-dependence of the higher twists h N (x) have been extracted from the data in a 
model independent way. It is shown that the size of their contribution to g\ is not 
negligible and their shape depends on the target: h p (x) changes sign while h n (x) is a 
non-negative function in the measured x region. 

We have found that the polarized parton densities depend on whether the higher 
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twist terms are or are not included in the analysis of g\. Moreover, the NLO polarized 
parton densities extracted from the g\ data in the presence of higher twist terms are 
in good agreement with those determined by our previous fits [5] to the data on g\jF\ 
and A\ using for the structure functions g\ and F\ only their leading twist expressions 
in NLO QCD. This observation confirms once more that the higher twist corrections 
to gi/F 1 and A 1 are negligible so that in the analysis of g\jF x and A x data it is enough 
to account only for the leading twist of the structure functions g\ and Fi. On the 
other hand, in fits to the g\ data themselves the higher twist contribution to g± must 
be taken into account. The latter is especially important for the LO QCD analysis of 
the inclusive and semi-inclusive DIS data. 
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Table 1. Parameters of the LO and NLO(JET) input parton densities at 
Q 2 = 1 GeV 2 as obtained from the best fits to the world g\/Fi and data 
including the HT corrections to g\. The errors shown are total (statistical and 
systematic). The parameters marked by (*) are fixed. Note that the TMC are 
included in ((?i)lt- 



Fit 


(9i)lo + h(x)/Q 2 


(9i)nlo + h(x)/Q 2 


DF 


185 - 16 


185 - 16 


x 2 


150.7 




145.0 


X 2 /DF 


0.892 




0.858 




0.926* 




0.926* 


a u 


0.000 ± 0.002 


0.312 ± 0.048 


In 


1.556 ± 0.261 




0* 


Vd 


- 0.341* 




-0.341* 


ad 


0.000 ± 0.005 


0.000 ± 0.049 


Id 


2.808 ± 1.249 




0* 


Vs 


- 0.072 ± 0.008 


- 0.045 ± 0.007 


a s 


0.601 ± 0.064 


1.5 


83 ± 0.434 


% 


0.803* 


0.8 


03 ± 0.244 


a g 


0.376* 


0.376 ± 0.503 


Xi 




) [GeV 2 ] 




0.028 


0.013 ± 0.036 


0.064 ± 0.044 


0.100 


- 0.076 ± 0.032 


- 0.007 ± 0.034 


0.200 


- 0.145 ± 0.032 


- 0.060 ± 0.035 


0.350 


- 0.030 ± 0.035 


- 0.008 ± 0.038 


0.600 


0.035 ± 0.019 


0.026 ± 0.021 


Xi 


h n ( Xi 


) [GeV 2 ' 




0.028 


0.234 ± 0.073 


0.178 ± 0.078 


0.100 


0.192 ± 0.048 


0.199 ± 0.050 


0.200 


0.035 ± 0.056 


0.079 ± 0.059 


0.325 


0.072 ± 0.071 


0.055 ± 0.073 


0.500 


0.023 ± 0.043 


- 0.020 ± 0.040 
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Table 2. The values of x 2 for the LO and NLO QCD fits without HT included 
compared to those when the HT corrections to g 1 are taken into account. The 
TMC are included in (#i)lt • 



Fit 


LO(HT=0) 


NLO(HT=0) 


LO+HT 


NLO+HT 


x 2 


244.5 


218.8 


150.9 


145.0 


DF 


185 - 6 


185 - 6 


185 - 16 


185 - 16 


X 2 /DF 


1.36 


1.22 


0.893 


0.858 
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Figure Captions 



Fig.l. Effective higher twist contribution h Al (x) to the spin asymmetry A±(x, Q 2 ) 
extracted from the data. Compared to our 1999 result (Fig. 4 [3] and Fig. 1 [5]) 
we present here the results of a new analysis including in the data set the SLAC/155 
proton data not available at that time. 

Fig. 2. Higher twist corrections to the proton and neutron g\ structure functions 
extracted from the data on gi in the case of LO and NLO QCD approximation for 
9i(x,Q 2 ) LT - 

Fig. 3. Higher twist corrections to the proton and neutron g\ structure functions 
extracted from the data when the leading twist calculations of gi(x, Q 2 )nlo are per- 
formed in different factorization schemes. 

Fig. 4. NLO( JET) polarized parton densities PD(gf LO + HT) (solid curves) compared 
to PD(gf LO /F^ LO ) (dashed curves) at Q 2 = 1 GeV 2 (see the text). 

Fig. 5. LO polarized parton densities PD(gl° + HT) (solid curves) compared to 
PD(^} j0 /F 1 L0 ) (dashed curves) at Q 2 = 1 GeV 2 (see the text). 

Fig. 6. Comparison of the proton structure function g\° calculated using the polarized 
parton densities PD(^° + HT) (solid curve) and PD(g\° / F^ ) (dashed curve) with 
SLAC/E143 proton data. 
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